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Abstract 

The three-dimensional instability of two coupled electromagnetic waves in an unmagnetized 
plasma is investigated theoretically and numerically. In the regime of two-plasmon decay, where 
one pump wave frequency is approximately twice the electron plasma frequency, we find that the 
coupled pump waves give rise to enhanced instability with wave vectors between those of the 
two beams. In the case of ion parametric decay instability, where the pump wave decays into 
one Langmuir wave and one ion acoustic wave, the instability regions are added with no distinct 
amplification. Our investigation can be useful in interpreting laser-plasma as well as ionospheric 
heating experiments. 
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I. INTRODUCTION 



The development of high-power lasers has made it possible to study the relativistic regime 
of laser-plasma interactions. This is of fundamental importance for laser-induced heating 
and inertial confinement fusion [1] and poses also the possibility to test fundamental physics. 
The first theoretical studies of relativistic waves in plasmas date back to the fifties [2J, where 
it was recognized that the the quivering velocity of electrons may lead to relativistic mass 
increase in an ultra-intense electromagnetic wave. The relativistic effects can result in self- 
modulation and self-focusing of electromagnetic waves in plasmas [3] . The instability of rel- 
ativistically large amplitude electromagnetic waves has been studied for magnetized electron 
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The instability of relativistically strong 



laser light in an unmagnetized plasma has then been revisited . The presence of mu 



laser beams in a plasma can give rise to a new set of interesting phenomena 



3,Q, 



14 



tip 



15| 



An important application of two interacting laser beams in plasmas is the excitation of large 



amplitude Langmuir waves, which in turn accelerate electrons to ultra-relativistic speeds [13]. 
The nonlinear coupling between two electromagnetic waves in plasmas can be described by 
a system of coupled nonlinear Schrodinger equations that model nonlinear interactions be- 



tween localized light [15|, 



161 ] and Langmuir or ion-acoustic waves. At strongly relativistic 



intensities, laser beams can give rise to fast plasma waves via higher-order nonlinearities 



12 



13 



171], or via the beat wave excitation at 



requencies different from the electron plasma 
frequency 18] . Particle-in-cell simulations [19j show that large- amplitude electron plasma 
waves can be excited by colliding laser pulses, or by two co-propagating electromagnetic 
pulses where a long trailing pulse is modulated efficiently by the periodic plasma wake be- 
hind the heading short laser pulse [201 ] . The effects on parametric instabilities of a partially 
incoherent pump wave (with a distribution of wave modes) was investigated both theoret- 
ically 21] and experimentally [221] , where it was found that the effect of finite bandwidth 
is, in general, to increase the instability thresholds and lower the growth rate. The nonlin- 



ear interaction 
considered 23, 



aetween two electromagnetic waves in the Earth's ionosphere has also been 
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251 ] . In that case, it is, in addition, important to focus attention on the 



collisional coupling between the waves, e.g. 



semiconductor plasmas [29| 
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281 ] . Related phenomena occur also in 
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Our previous treatment of the instability of two coupled laser beams in an electron-ion 



plasma 



T 

3JJ and a two-temperature electron plasma |32| were limited to the investigation 



of the relativistic Raman and Brillouin scattering instabilities in one dimension, and two- 
dimensional effects were only partly included. The purpose of the present paper is to include 
the multi-dimensional effects, which are particularly important for the cases where an elec- 
tromagnetic wave decays into one low-frequency wave and one high-frequency electrostatic 
wave that are propagating obliquely to the pump wave. 



II. MODEL EQUATIONS 



Let us consider the propagation of intense laser light in an electron-ion plasma. The 
dynamics of the high frequency laser light is governed by 
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|Aft| 2 )Aft = 0, (1) 



dt 2 

where uj pe = (Aimoe 2 /m e ) 1//2 is the electron plasma frequency, n is the equilibrium number 
density, e is the magnitude of the electron charge, m e is the electron mass, v x e is the electron 
thermal speed, c is the speed of light, and the normalized slow time-scale electron number 
density perturbation is N es = n es /n . The latter is excited by the ponderomotive force of 



the high frequency waves. If the ions are considered as immobile, we have 
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whereas if the electrons are treated as inertialess 31] 

V ot J m e mjC A 
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where c s is the sound speed. 

We will here consider two large amplitude electromagnetic waves Aft = Ai + A2. We 
then have |A h | 2 = |Ai| 2 + |A 2 | 2 + 2Ai • A 2 . 

Using equations (l)-(3) we can now derive a nonlinear dispersion relation. The calcu- 
lations are straightforward but lengthy. Following closely the analysis of Ref. j^ij], it then 
turns out to be convenient to introduce a characteristic velocity v t s that is defined by 
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where u>j± = to ± u>oj, kj± = k ± k j, uoj and k j is the pump frequency and wavevector 
of the pump wave j (j = 1, 2), whereas lo and k are associated with the electrostatic 
low-frequency fluctuations. Furthermore, we have here introduced the pump wave quiver 



velocity v 0j - = eAqJm e c, Dj± 
sideband damping [33j, and &j±£j± 
use wgj- = + fcgj-C 2 . 
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— k 2 ±c 2 + iuij±jj± where jj± represents the 
v pe ~ ^j± v Te- For the pump frequencies, we 



In the present paper we will, for simplicity, further assume that Ai ■ A2 ~ 0. This 



means that double resonance parametric phenomena [34| . where the difference between the 
frequencies cg> i and CU02 is close to twice a natural frequency, will be neglected. Choosing 
two transverse waves that propagate in the y— and z-directions, respectively with the pump 
velocities in the z— and y-directions, respectively, we note that (4) reduces to 
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With these limitations, and following Refs. 31] and 33j, the nonlinear dispersion relation 
turns out to be 
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In Eq. (6), as well as below, we have normalized A j by m e c 2 /e. For electron Langmuir 
waves, we have [3JJ 

Ql = < i 1 ~ u 2 -3k 2 v 2 -^ ■ (?) 



pe , 



whereas for ion acoustic waves [31] 

Qia = 
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III. RESULTS 



We have carried out a numerical analysis of Eq. ([6]), where we have assumed that the 
frequency is complex valued, where the imaginary part of lo represents the growth rate. In 
our treatment, we have concentrated on three- wave decay processes where the important 
terms in Eq. are the ones with down-shifted daughter waves. Hence, we have kept the 
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FIG. 1: The growth rate (normalized by oj pe ) of the two-plasmon decay of one laser beam for 
different electron thermal speeds a) vx e = 0, b) vx e = 0.005 c, c) vx e = 0.0075 c, and d) vx e = 0.01 c. 
The pump amplitude is Aoi = 0.01 z, the pump wavevector koi = \/3yujp e /c. The second laser 
beam intensity is set to zero (A02 = 0). 




FIG. 2: The growth rate (normalized by uj pe ) of the two-plasmon decay of one laser beam for 
wavevectors a) koi = 1.65yw pe /c, b) koi = V3yaj pe /c, c) koi = 2yw pe /c, and d) koi = 2.5ycj pe /c. 
The pump amplitude is Aoi = 0.01 z, and the electron thermal speed VTe = 0.01c. The second 
laser beam intensity is set to zero (A02 = 0). 
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FIG. 3: The growth rate (normalized by uj pe ) of the two-plasmon decay of two coupled laser 
beams for wavevectors a) koi = 1.65yu pe /c and ko2 = 1.65zuj pe /c, b) koi = ^/3yu) pe /c and ko2 = 
\/3zLijp e /c, c) koi = 2yuj pe /c and ko2 = 2zoj pe /c, and d) koi = 2yio pe jz and ko2 = \/3zw pe /c. 
The pump amplitudes are Aoi = 0.01 z and A02 = 0.01 y, and the electron thermal speed is 

VTe = 0.01 C. 

resonant terms with subscripts "-" in Eq. (JBj) but neglected those with subscripts "+". For 
large wavenumbers, we recover almost identically the stimulated Raman and Brillouin cases 
treated in Ref. 3JJ. Thus, we shall here concentrate on the regimes of smaller wavenumbers 
in which the two-plasmon decay and ion parametric decay instabilities become important, 
and where it is crucial to have a fully multi-dimensional treatment. 

We first consider the two-plasmon case in which the low-frequency wave is a Langmuir 
wave, and where Q in Eq. (jSJ) equals Ql in ([Z|). The growth rates (normalized by u pe ) are 
presented in Figs. 1-3. We have here denoted the unit vectors in the y and z directions 
by y and z. In Fig. 1, the electron thermal effects are investigated. The dispersion of the 
Langmuir waves due to the electron pressure causes the instability to restricted in a bounded 
domain in wavevector space, with a maximum growth rate in a direction perpendicular to the 
propagation direction of the laser beam. In Fig. 2, we have investigated how the instability 
depends on the pump wavenumber of a single electromagnetic beam. We find that the 
instability is dominant for wavenumbers equal to or larger than \/3u pe /c, and that the 
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FIG. 4: The growth rate (normalized by ui pe ) of the ion parametric decay instability of one laser 
beam for different electron thermal speeds a) VTe = 0, b) vx e = 0.005 c, c) vx e = 0.0075 c, and d) 
VTe = 0.01 c; in each case the ion acoustic speed is set to c s = 0.006 vx e - The pump amplitude is 
Aoi = 0.01 z, and the pump wavevector is koi = 0.1yw pe /c. The second laser beam intensity is set 
to zero (A02 = 0). 

instability vanishes for smaller wavenumbers of the pump wave. The maximum instability 
occurs for wavevectors almost perpendicular to the laser beam propagation direction. The 
instability maximum occurs for wavenumbers much larger than the pump wavenumber. 
Eventually kinetic effects will become important and electron Landau damping will decrease 
the growth rate. The case of two coupled electromagnetic beams is investigated in Fig. 3. 
While beam 1 is directed in the y-direction as in Figs. 1 and 2, beam 2 is here in the 
z-direction, perpendicular to beam 1. We see in panel a) of Fig. 3 that the coupled beam 
system gives rise to a new instability with a maximum growth rate in the direction of the 
dichotome in the center between the two beam propagation directions. For wavenumbers 
larger than ^u pe /c, the instability becomes more evenly distributed in all directions, but 
with a well-defined maximum growth rate for some wavenumber. For the case where the 
wavenumber of beam 2 is smaller than the one of beam 1, in panel d), we see a superposition 
of the two instability regions for the separate beams. We note that similar effects can also 
appear due to direct subharmonic wave generation in nonuniform plasmas 35]. 
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FIG. 5: The growth rate (normalized by ui pe ) of the ion parametric decay instability of one laser 
beam for wavevectors a) koi = 0, b) koi = 0.1yuj pe /c, c) koi = 0.5y<^ pe /c, and d) koi = yu pe /c. 
The pump amplitude is Aoi = 0.01 z, the electron thermal speed is vt& = 0.01c, and the ion 
acoustic speed is c s = 6 x 10 _5 c. The second laser beam intensity is set to zero (A02 = 0). 

We next investigate the parametric decay instability in which the electromagnetic wave 
decays into one electrostatic wave and one low-frequency ion acoustic wave, where Q in 
Eq. (JHJ) equals Qia in (EJ). The growth rates are presented in Figs. (4)-(6). We have here 
concentrated on the long wavelength limit where the main instability is the ion parametric 
decay instability, in which the electromagnetic wave decays into one slightly frequency- 
downshifted electrostatic wave and one low-frequency ion acoustic wave. In Fig. 4, we have 
studied the thermal effect on the instability of one single electromagnetic beam. We see 
that for higher electron thermal and ion acoustic speeds, the region of instability becomes 
smaller in wavenumber space, and that there is a well-defined maximum of the instability 
for propagation almost perpendicular to the wave propagation direction. In Fig. 5, we have 
considered different wavenumbers of the pump wave. We see that for large wavenumbers, the 
maximum instability occurs for larger wavevectors perpendicular to the pump wavevector. 
Finally, we consider the interaction between two coupled electromagnetic beams in Fig. 6. 
For the cases of equal pump amplitudes and lengths of the wavevectors, in panels a)— c), the 
instability becomes almost rotationally symmetric, with equal maximum growth rates in all 
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FIG. 6: The growth rate (normalized by to pe ) of the ion parametric decay instability of two coupled 
laser beams for wavevectors a) koi = 0.1yu;p e /c and ko2 = 0.1zw pe /c, b) koi = 0.5yw pe /c and 
ko2 = 0.5zw pe /c, c) koi = yujp e /c and ko2 = zw pe /c, and d) koi = yu) pe /c and ko2 = 0.5zw pe /c. 
The pump amplitudes are Aoi = 0.01 z and A02 = 0.01 y, the electron thermal speed is vr e = 0.01 c 
and the ion acoustic speed is c s = 6 x 10~ 5 c. 

directions. For the case where the wavenumber of beam 2 is smaller than the one of beam 
1, in panel d), we see a superposition of the two instability regions for the separate beams. 
For the ion parametric decay instability, we could not see the same distinct amplification of 
the instabilities as we could see for the two-plasmon decay in panel a) of Fig. 3. 

In conclusion, we have investigated the instability of two coupled large-amplitude elec- 
tromagnetic waves in a plasma. Our investigation shows that two-plasmon decay plays an 
important role for pump frequencies approximately two times larger than the background 
electron plasma frequency, and that one has a maximum growth rate perpendicular to the 
propagation direction of the electromagnetic wave. For two coupled electromagnetic beams, 
there is here a new and stronger instability in a direction between the two laser beams. 
The ion parametric decay instability, in which an electromagnetic wave decays into one 
Langmuir wave and one ion acoustic wave, is important for long wavelengths of the electro- 
magnetic pump wave, where the instability leads to ion acoustic waves that are propagating 
perpendicularly to the electromagnetic beam direction. Here, the addition of a second elec- 
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tromagnetic beam leads to a superposition of the instabilities. However, we see no strong 
nonlinear amplification due to the two beams. Our study could be important for both 
laser-plasma interactions and for ionospheric heating experiments [36l. 1371. l38l. 139]. 
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